
BELLE 



Belle Preprint 2002-28 
KEK Preprint 2002-83 
TIT-HPE-02-01 



(N 
> 



Measurement of the Electroweak Penguin Process 

(N"' B Xs£+£- 

o : 

^ : J. Kaneko,^^ K. Abe,^ K. Abe,^^ T. Abe,^^ I. Adachi,^ Byoung Sup Ahn/^ H. Aihara,^^ 

M. Akatsu,23 Y. Asano,^^ T. Aso,^^ V. Aulchenko,^ T. Aushev/^ A. M. Bakich,^^ Y. Ban,^^ 
E. Banas,28 W. Bartel,^ A. Bay/^ P. K. Behera,^^ a. Bondar,^ A. Bozek,^^ M. Bracko,^!'^^ 
J. Brodzicka,28 T. E. Browder,^ B. C. K. Casey,^ P. Chang,^^ Y. Chao,^^ K.-F. Chen,^^ 
^ ; B. G. Cheon,^o R. Chistov/^ Y. Choi,^° Y. K. Choi,^° M. Danilov/^ L. Y. Dong," 

S. Eidelman,^ V. Eiges/^ Y. Enari,^^ C. W. Everton,^^ F. Fang,^ H. Fujii,^ C. Fukunaga,^^ 
N. Gabyshev,9 A. Garmash,^'^ T. Gershon,^ R. Guo,^^ J. Haba,^ F. Handa,^^ 
0\ '. T. Hara,^2 y. Harada,^^ N. C. Hastings,^^ H. Hayashii,^^ M. Hazumi,^ E. M. Heenan,^^ 

g : 1. Higuchi,^^ T. Higuchi,^^ L. Hinz/^ T. Hojo,32 Y. Hoshi,^^ W.-S. Hou,^^ H.-C. Huang,^^ 

OO ! T. Igaki,^^ Y. Igarashi,^ T. lijima,^^ K. Inami,^^ A. Ishikawa,^^ H. Ishino,^^ R. Itoh,^ 

^ ; H. Iwasaki,^ Y. Iwasaki,^ H. K. Jang,^^ J. H. Kang,^^ J. S. Kang/^ N. Katayama,^ 

O ' H. Kawai,^ Y. Kawakami,^^ N. Kawamura/ T. Kawasaki, H. Kichimi,^ D. W. Kim,^° 

■ Heejong Kim,^^ H. J. Kim,^^ H. O. Kim,'^° Hyunwoo Kim,^^ S. K. Kim,^^ K. Kinoshita,^ 

^ • S. Kobayashi,^^ S. Korpar,^^'^^ P. Krizan,^°' P. Krokovny,^ R. Kulasiri,^ S. Kumar, 

gn! Y.-J. Kwon,55 J. S. Lange,^'^^ G. Leder,^^ g ^ 39 j 38 ^ Limosani,^^ R.-S. Lu,^^ 

^ : J. MacNaughton,^2 q Majumder,*^^ p_ Mand^^ jy Marlow,^^ S. Matsumoto,"^ 

> : T. Matsumoto,^^ W. Mitaroff,^^ k. Miyabayashi,^^ Y. Miyabayashi,^^ H. Miyake,^^ 

'O : G. R. Moloney,22 T. Mori,^ T. Nagamine,^^ Y. Nagasaka,^° T. Nakadaira,^^ E. Nakano,^^ 

^ ; M. Nakao,9 J. W. Nam,^^ k. Neichi,^^ S. Nishida,i^ O. Nitoh,^^ S. Noguchi,^^ T. Nozaki,^ 

^ ' S. Ogawa,^3 rj Ohshima,23 T. Okabe,^^ S. Okuno,!^ S. L. Olsen,^ Y. Onuki,^^ 

W. Ostrowicz,28 H. Ozaki,^ P. Pakhlov,^^ ^ Palka,^^ C. W. Park,!^ H. Park,!^ K. S. Park,^o 
J.-P. Perroud,^^ M. Peters,^ L. E. Piilonen,^^ Root,^ K. Rybicki,^^ H. Sagawa,^ 
Y. Sakai,^ H. Sakamoto,^'' M. Satapathy,^^ A. Satpathy,^'^ O. Schneider, S. Schrenk,^ 
C. Schwanda,^' S. Semenov,^^ K. Senyo,^^ R. Seuster,^ H. Shibuya,"^^ B. Shwartz,^ 
V. Sidorov,2 J. B. Singh,^^ N. Soni,33 S. Stanic,^^'!] K. Sumisawa,^ T. Sumiyoshi,^^ 
K. Suzuki,^ S. Suzuki, S. Y. Suzuki,^ S. K. Swain, ^ H. Tajima,''^ T. Takahashi,^^ 
F. Takasaki,^ K. Tamai,^ N. Tamura,^° J. Tanaka,^^ M. Tanaka,^ G. N. Taylor,^^ 
Y. Teramoto,^^ S. Tokuda,^^ M. Tomoto,^ T. Tomura,^^ K. Trabelsi,^ W. Trischuk,^^'! 

T. Tsuboyama,^ T. Tsukamoto,^ S. Uehara,^ K. Ueno,^^ S. Uno,^ Y. Ushiroda,^ 
G. Varner,^ K. E. Varvell,^^ C. C. Wang,^^ C. H. Wang,^^ J. G. Wang,^^ ^.-Z. Wang,^^ 
Y. Watanabe,'''' E. Won,^^ B. D. Yabsley,^^ Y. Yamada,^ A. Yamaguchi,^^ Y. Yamashita,^^ 
M. Yamauchi,^ H. Yanai,^° J. Yashima,^ M. Yokoyama,^^ Y. Yuan,^^ Y. Yusa,'*^ 
C. C. Zhang," J. Zhang,^^ Z. P. Zhang,^^ Y. Zheng,^ V. Zhihch,^ and D. Zontar^^ 



1 



(The Belle Collaboration) 

^ Aomori University, Aomori 
^Budker Institute of Nuclear Physics, Novosibirsk 
^Chiha University, Chiba 
^Chuo University, Tokyo 
^University of Cincinnati, Cincinnati OH 
^Deutsches Elektronen-Synchrotron, Hamburg 
University of Frankfurt, Frankfurt 
^University of Hawaii, Honolulu HI 
^High Energy Accelerator Research Organization (KEK), Tsukuba 
Hiroshima Institute of Technology, Hiroshima 
Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 
^^Institute of High Energy Physics, Vienna 
^^Institute for Theoretical and Experimental Physics, Moscow 
^'^J. Stefan Institute, Ljubljana 
^^Kanagawa University, Yokohama 
^^Korea University, Seoul 
^'^Kyoto University, Kyoto 
^^Kyungpook National University, Taegu 
^^Institut de Physique des Hautes Energies, Universite de Lausanne, Lausanne 

University of Ljubljana, Ljubljana 

University of Maribor, Maribor 
University of Melbourne, Victoria 
^^Nagoya University, Nagoya 
^^Nara Women's University, Nara 
National Kaohsiung Normal University, Kaohsiung 
^^National Lien-Ho Institute of Technology, Miao Li 
^'^ National Taiwan University, Taipei 
^^H. Niewodniczanski Institute of Nuclear Physics, Krakow 
^^Nihon Dental College, Niigata 
^^Niigata University, Niigata 
Osaka City University, Osaka 
Osaka University, Osaka 
^^Panjab University, Chandigarh 

^^Peking University, Beijing 
Princeton University, Princeton NJ 
^^RIKEN BNL Research Center, Brookhaven NY 

^'^Saga University, Saga 
University of Science and Technology of China, Hefei 
^^Seoul National University, Seoul 
Sungkyunkwan University, Suwon 
'^^ University of Sydney, Sydney NSW 
^^Tata Institute of Fundamental Research, Bombay 
^^Toho University, Funabashi 
^^Tohoku Cakuin University, Tagajo 
^^Tohoku University, Sendai 



2 



University of Tokyo, Tokyo 
'^"'Tokyo Institute of Technology, Tokyo 
Tokyo Metropolitan University, Tokyo 
Tokyo University of Agriculture and Technology, Tokyo 
^'^Toyama National College of Maritime Technology, Toyama 
University of Tsukuba, Tsukuba 
^^Utkal University, Bhubaneswer 
Virginia Polytechnic Institute and State University, Blacksburg VA 
"'^ Yokkaichi University, Yokkaichi 
^^Yonsei University, Seoul 
(Dated: August 18, 2002) 

Abstract 

We report the first measurement of tlie brandling fraction for the inclusive decay B Xsi'^i~ , 
where i is either an electron or a muon and Xg is a hadronic recoil system that contains an s-quark. 
We analyzed a data sample of 65.4 x 10^ B meson pairs collected with the Belle detector at the 
KEKB e+e- asymmetric-energy collider. We find B{B Xsi+£~) = (6.1 ± 1.4(stat) li:^(syst)) x 
10~^ for dilepton masses greater than 0.2 GeV/c^. 

PACS numbers: 13.20.He, 14.65.Fy, 14.40.Nd 
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The first observation of tlie flavor-clianging-neutral-current (FCNC) weak decay process 
B — >■ Ki~^i~, recently reported by Belle [|l|, opens a new window for searches for physics 
beyond the Standard Model (SM) [Q]. There are no SM first-order weak decay processes that 
can produce such decays; the dominant SM processes are second-order 6-quark to s-quark 
processes with a virtual t-quark in a loop (electroweak penguin) or a box diagram. Non-SM 
processes could produce sizable modifications to the decay amplitude due to contributions 
from virtual non-SM particles (such as charged Higgs or SUSY particles) in the loop. 

The decay amplitude is often described by an effective Hamiltonian in which non-SM 
contributions can modify the Wilson coefficients C7, Cg and Cio- The magnitude of allowable 
modifications of C7 is constrained by the measured rate for inclusive B —>■ Xgj decays 
IP, ^. The Belle measurement of the exclusive process B Ki^i~ has been used to 
determine the best limits on non-SM contributions to Cg and Ciq 0. The usefulness of 
exclusive measurements is limited by the large theoretical uncertainties associated with the 
s-quark to K meson hadronization process. Since these uncertainties are not as severe for 
inclusive processes, measurements of the branching fraction and the dilepton and hadronic 
mass spectra for B —>■ Xsi~^i~ will provide more stringent and less model-dependent probes 
for new physics. This process has not been measured; CLEO has reported a 90% confidence 
upper limit of B{B ^ XJ+i') < 4.2 x IQ-^ 0. 

In this Letter, we present the results of a measurement of the branching fraction for the 
inclusive decay B — >■ Xsi~^i~, where B is either B^ or B^, i is either an electron or a muon, 
and Xs is a hadronic recoil system that contains an s-quark. Here and throughout the paper, 
the inclusion of charge conjugated modes is implied. We use a data sample collected at the 
T(45') resonance with the Belle detector at the KEKB e'^e~ asymmetric-energy collider (3.5 
GeV on 8 GeV) ||^. This sample contains (65.4 ± 0.5) x 10^ B meson pairs, corresponding 
to an integrated luminosity of 60 fb~^. 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50- layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), time-of-flight scintillation counters (TOF), and 
an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented with resistive plate counters to identify muons (KLM). 
The detector is described in detail elsewhere p. 

We reconstruct charged particle trajectories with the CDC and SVD. Electron identifica- 
tion is based on the position and shower shape of the cluster in the ECL, ratio of the cluster 
energy to the track momentum (E/p), specific energy-loss measurement (dE/dx) with the 
CDC and the response from the ACC. Muon identification is based on the hit positions 
and the depth of penetration into the ECL and KLM. We require the electron's (muon's) 
laboratory momentum to be greater than 0.5 GeV/c (1 GeV/c). We find an electron (muon) 
selection efficiency of 92.5% (91.3%) with a (0.2 ± 0.06)% ((1.4 ± 0.04)%) pion to electron 
(muon) mis-identification probability. Charged kaon candidates are selected by using in- 
formation from the ACC, TOF and dE/dx in the CDC. The kaon selection efficiency is 
90% with a pion to kaon mis- identification probability of 6%. The remaining charged tracks 
are assumed to be pions. We select — >■ tt^tt" candidates with invariant mass within 
15 MeV/c^ of the mass. We impose additional Kg selection criteria based on the dis- 
tance and the direction of the Kg vertex and the impact parameters of daughter tracks. We 
require the charged tracks other than those used in the Kg reconstruction to have impact 
parameters with respect to the nominal interaction point of less than 0.5 cm in the radial 
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direction and 3 cm along the beam direction. 

We reconstruct photons from ECL energy clusters that have no associated charged tracks. 
The shower shape is required to be consistent with an electromagnetic cluster and the energy 
to be greater than 50 MeV. We reconstruct vr^ — > 77 candidates from photon pairs with 
invariant mass within 10 MeV/c^ of the nominal 7r° mass. 

The Xg system is reconstructed in 18 different combinations of either a or Kg com- 
bined with to 4 pions, of which up to one 7r° is allowed. We combine the Xg with two 
oppositely charged leptons to form a B candidate. We identify the B — > Xsi~^i~ signal 



with the beam-energy constrained mass, Mbc = \J {E^^^/ c^Y — (p^^/c)^, where -Ebeam 
p'^ are the beam energy and the B candidate momentum calculated in the center-of-mass 
(CM) system, respectively. We find the average Mbc resolution is dbc = 2.8 MeV/c^. The 
energy difference /S.E = E'^^ — -Ebeam' where E'^^ is the CM energy of the B candidate, is 
combined with other variables to provide background suppression. 

To optimize the selection criteria and to determine their signal efficiencies, we use the 
following signal model. The SM calculation of Ref. is used for the dilepton mass {M^+i-) 
spectrum. We require Mi+£- to be greater than 0.2 GeV/c^ to remove the virtual photon 
contribution, b — > 57* — > se'^e^ |jT^. We model the recoil mass {MxJ spectrum using a non- 
resonant Fermi motion model [rT| and the JETSET hadronization program |T^. For Mx^ < 
1.1 GeV/c^, however, the spectrum is modeled by the sum of the exclusive B — > K^*H'^i~ 
components as predicted by the SM |jl3|. Using this model we find that our reconstructed 



final states account for 81% of the total Xg states. Figures |l|(a) and (b) show the M£+£- and 
spectra from this model. 
There are two background sources that can peak in Mbc and AE. The first is hadronic B 
decays into one kaon plus multiple pions {B XsTt+tt^) from abundant decays such as -B — *■ 
D^*^mi (n > 1). We estimate this background contribution by reconstructing B — > X^Tr'^Tr" 
events without the lepton identification criteria and multiplying by the measured momentum 
dependent pion mis-identification probability. We find contributions of 2.6 ± 0.2 events to 
the B XsyU+yU^ signal and 0.1 ± 0.05 events to B ^ Xse^e~ . We refer to this as the 
fake background, and subtract it from the signal yield. The second is from B ^ J / ipXg and 
B — > ip'Xg where J/ip and ip' decay into dileptons. These decay modes have the same final 
states and, in principle, interfere with the signal. For this study, these charmonium decays 
are explicitly vetoed. The veto windows are —0.6 to +0.2 GeV/c^ (—0.35 to +0.2 GeV/c^) 
around the J/ip mass for the e^e~ (/i^/i~) channel, and —0.3 to +0.15 GeV/c^ around the 
ijj' mass for both channels. We find the background from this source is negligibly small using 
a 5 — *• J/ipXs Monte Carlo (MC) sample in which the Mn+n- spectrum reproduces that of 
data. 

The largest background sources are random combinations of dileptons with a kaon and 
pions that originate from continuum qq [q = u, d, s, c) production or from semileptonic B 
decays. We reject 83% of the qq background with a signal efficiency of 90% by using a 
Fisher discriminant [1^ (^cont) based on a modified set of Fox- Wolfram moments [15| that 
differentiate the event topology. In the semileptonic B decay background, both B mesons 
decay into leptons or two leptons are produced from the b c —>■ s,d decay chain. We 



combine the missing mass (Mmiss) and the total visible energy (-Evis) [16| into another Fisher 
discriminant (^si) to reject 85% of events with two neutrinos with a signal efficiency of 91%. 

We further reduce the backgrounds using AE and the cosine of the B flight direction 
(cos^b) with respect to the e~ beam direction in the CM frame. First, we select events 
with \AE\ < 40 MeV (~3(Ta£;). We then calculate likelihoods Cs,b = Psb ^ pTb^ where 
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j9^B and p's'^B uncorrelated probability density functions (PDF) for AE and cos 9b 

for the signal (S) and the background (B), respectively, and form a likelihood ratio CTZ = 
Csl{Cs + Cb)- We model with a Gaussian primarily determined from a signal MC 
sample and calibrated using B J/ifjXa data; is modeled with a linear function with 
a slope determined from a MC sample that contains 6 — c decays and gg events. The signal 
follows a 1 — cos"^ 6b distribution, while the background distributes uniformly in cosOb- By 
requiring CTZ > 0.8, we retain 75% of the signal while removing 80% of backgrounds. 

For events with multiple candidates that pass the selection criteria, we choose the combi- 
nation with the largest value of CTZ. We find that the correct candidate is reconstructed in 
73% of the events. We reject candidates with Xg invariant mass greater than 2.1 GeV/c^. 
This condition removes a large fraction of combinatorial background while retaining (93±5)% 
of the signal. 

Signal MC samples are used to determine reconstruction efficiencies of (3.7 ±0.4 ±0.5)%, 
where the first error is systematic and the second error is due to the model uncertainty. 
Here, an equal production rate is assumed for B^B° and B^B~ . Systematic errors include 
uncertainties from the tracking efficiency (2.0% per track), reconstruction efficiency 
(8.7% per K^), vr" reconstruction efficiency (6.8% per vr"), electron identification (1.8% 
per electron), muon identification (2.2% per muon), identification (2.5% per K^)^ 7r+ 
identification (0.8% per tt"*"), and the background suppression criteria (3.0% for J^cont, 
and CTZ, combined, estimated using B J/ipXg data). 

The model uncertainty includes the following sources. The largest error is due to un- 
certainties in the SM branching fraction of the exclusive modes (11%). The uncertainty 
in the Mx^ spectrum (4%) is due to the Fermi momentum (pp) and spectator quark mass 
{rriq) parameters. We take pp = 0.4 GeV/c and niq = GeV/c^, and vary the parameters 
over a range allowed by the measured heavy quark effective theory (HQET) parameters Ai 
and A [1^. The uncertainty in the fraction of the unmeasured modes (2.1%) and the 
uncertainty due to the fractions of 7r° and Kg contained in Xg (4.6%) are estimated by 
comparing the inclusive hadron (vr^. Kg, rj, 0, etc.) production rates between continuum 
data and JETSET. 

We determine the signal yield from an unbinned maximum likelihood fit to the 
distribution as shown in Fig. 1^. We model the signal with a Gaussian function and the 



background with a threshold function |[T8[. The width of the signal Gaussian is obtained 
from B — »• J/ijjXs data. The background shape is obtained from the background MC 
sample. We verify that the threshold function obtained from MC reproduces the shape of 
data for B X^e/x combinations (Fig. 0(d)). The fit results are given in Table |; we find 
60.1 ± 13.9(stat)l5:4(syst) events for the combined B Xgi+i" modes (Fig. |(c)) with a 
statistical significance of 5.4(7. Here, the systematic error is obtained from the maximum 
deviation when the background and signal shapes are varied by one standard deviation of the 



statistical error in their shape parameters. The significance is defined as y — 2 ln(£o/'Cmax), 
where £max is the maximum likelihood and Co is the maximum likelihood when the signal 
yield is constrained to be zero. The fake background and the uncertainty in the background 
shape are included in the significance calculation. We calculate the branching fraction for 

Mi+e~ > 0.2 GeV/c2 to be 

B{B X,£+r) = (6.1 ± 1.4(stat) +}:t(syst)) x 10"^ 

where the systematic errors in the yield, efficiency, the number of B meson pairs and the 
model errors are added in quadrature to give the total systematic error. This can be com- 
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pared to the SM expectation B{B Xsi+i') = (4.2 ± 0.7) x 10"^ fig. Table | summarizes 
the branching fractions for B — >■ X^e^e" and B — > Xsfi'^fi~ separately, together with the 
number of candidates, signal yields, fake background estimations, efficiencies and the statis- 
tical significances. 

The dilepton mass spectrum is measured by dividing the data into M^+i- bins. For each 
bin, the signal yield is extracted from a fit to the Mbc distribution and the fake background 
contribution is subtracted. The result is shown in Fig. ^(c). Similarly, the recoil mass 
spectrum is obtained by dividing the data into Mx^ bins and extracting the signal yield for 
each bin. The result is shown in Fig. [^(d). With the current statistics, the Mi+i- and 
spectra are in agreement with SM expectations. Branching fractions for each bin are given 
in Table |l[ 

In summary, we present the first measurement of the inclusive branching fraction for 
the electroweak penguin decay B —>■ Xs£~^£~. The results are in agreement with the SM 
expectations and can be used to constrain extensions of the SM. 
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TABLE L Fit results for the number of candidates, signal yields, fake backgrounds, reconstruction 
efficiencies, statistical significances and branching fractions {B). Candidates are counted in a ±3o"bc 
window in My^^- 

mode candidates signal yield fake background efficiency (%) significance B{xlO^^) 
B^XsB+e- 96 25.5 ±11.2+11 0.1 ±0.05 3.9 ± 0.4 ± 0.5 3.4 5.0 ±2.3li 
B^Xsfi+fi-- 92 37.3 ±9.7l|j 2.6 ± 0.2 3.6 ± 0.4 ± 0.5 4.7 7.9 ±2.ll^ 
B^Xsl+i' 188 60.1 ±13.9t|^ 2.7±0.2 3.7±0.4±0.5 5.4 6.1 ±1.4li 



TABLE IL Branching fractions {B) for each bin of Me+^~ and Mxs- The first and second errors 
are statistical and systematic, respectively. 
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[1.0, 1.6] 


26.8 ± 9.5 
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11.1 ±5.4 
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26.8 ± 13.8115 
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FIG. 1: SM expectations for the (a) dilepton and (b) recoil mass spectra; the observed (c) dilepton 
and (d) recoil mass spectra (circles). Inner and outer error bars indicate the statistical and total 
errors, respectively. The histograms in (c), (d) show the SM expectations after all the selections 
are applied; histograms are normalized to the expected branching fractions. The gaps in (c) are 
due to the J ftp, V"' vetoes. The dotted line in (d) indicates the Mx^ < 2.1 GeV/c^ requirement. 
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FIG. 2: Beam-energy constrained mass distributions for (a) B XgC'^e^ , (b) B X^/x^yU^, (c) 
B Xgi'^l" and (d) B —>■ XgC^fi^. The solid lines indicate the fit results and the dotted lines 
show the sum of the background components. 
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